Parameter fluctuations in coupled chaotic systems 
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We study the effect of parameter fluctuations on synchronization of a coupled chaotic system. 
The fluctuations to the parameter can be random or it can be a periodic modulation. For random 
fluctuations we introduce a new quantity, the rate of fluctuations, apart from the statistical features 
of the fluctuations. Fluctuation rate in our study refers to the number of random modifications to 
parameters occurring in unit time. With a periodic modulation, the fluctuation rates can be the 
frequency of the modulating term. It is found that with high fluctuation rates the synchronization is 
stable irrespective of the statistical or mathematical features of fluctuation. It was also found that 
the low fluctuation rates can destroy synchronization even with a small amplitude. We analytically 
explain the observed phenomenon using the dynamical equations and numerically verified with 
a coupled system of Rossler attractor as an example. We also numerically quantify the relation 
between synchronization error and fluctuation rates. 

PACS numbers: 05.45.Ac, 05.45.Pq, 05.45.Xt, 05.45.Vx 



I. INTRODUCTION 

Synchronization of coupled chaotic systems has gen- 
erated a lot of research activities over the last sev- 
eral years. Synchronized behavior has been studied 
extensively in physical, chemical and biological sys- 
tems. Different types of synchronization such as complete, 
generalized, lag and phase synchrony are described in 
literature. One of the methods by which the synchro- 
nization of chaotic systems is achieved is by coupling 
two identical systems, which may be unidirectional or 
bidirectional[J 0, S HE 113 IH III El • Synchronization 
in arrays of coupled laser systems has also been investi- 
gated under various coupling schemes pi TToL fTTl ll2L Il3| . 
Complete synchronization of identical chaotic systems is 
also of considerable interest because of its applications 
in secure communication [T^, Il3| . By identical systems 
we mean a set of systems whose parameters are exactly 
equal. It is found that the synchronization is not robust 
when there is a small but finite mismatch of the param- 
eters of the systems 0, Q, @- In fact the effect of phase 
mismatch and a finite constant frequency detuning in a 
bidirectionally coupled Duffing oscillators is to destroy 
the synchronization altogether [l4j. 

In the present paper we address a different issue that 
is relevant in many practical physical systems. We study 
the effect of random fluctuations in the parameters of the 
system on the synchronization properties. Such a study 
is relevant and important since in a real physical sys- 
tem the parameters can often fluctuate randomly either 
due to some internal instabilities or due to some external 
perturbations. 



The paper is organized as follows; in section II. we 
consider random perturbations to one of the parameter 
that characterizes the synchronization. We present the 
criteria for the robustness of synchronization. In sec- 
tion III. the numerical results on two coupled Rossler 
systems with randomly fluctuating parameters are pre- 
sented. Section IV contains the discussions of the result. 



II. PARAMETER FLUCTUATIONS IN 
COUPLED SYSTEMS 

In this section we consider two identical dynamical sys- 
tems which are coupled together. To study the effect 
of fluctuations it is essential to identify one parameter 
whose mismatch is most effective in destroying synchro- 
nization. We denote this parameter as p. Then the equa- 
tions for the coupled systems are given by. 



X x = h( P i,X 1 )+Cf(X 2 -X 1 ) 
X 2 = fi(p2,X 2 )+Cf(X 1 -X 2 ) 



(1) 



Here C is the coupling constant. In reality it is difficult 
if not impossible to construct identical systems except in 
numerical simulations. This can also be due to the fact 
that the parameters could be fluctuating in time with a 
time scale of their own. We incorporate this by writing 
the parameter as 



Pi = Pq + £it 
P2=P0+ £,2t , 



(2) 
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where, £k and £2* are two delta correlated zero mean 
Gaussian random variables. A measure of the amplitude 
of fluctuations, we define Ap, as 



Ap = (| 6 P (t) |) t , 



(3) 



2 



where, Sp(t) = p\(t) — p 2 (t) and (■■■)t denotes time aver- 
age. 

To study the effect of time scales of parame- 
ter fluctuation, we define the fluctuation rate <f> — 
number of perturbuations / unit time. Different fluctua- 
tion rates can be achieved numerically by modifying the 
parameter as in Eqr|21only in certain chosen time steps. 
Rest of the time the value of the parameter remains con- 
stant at the modified value. The Error in synchrony is 
studied varying <fi. The effect of time scales has not been 
studied in literature and our results indicate that it is 
highly significant in determining the quality of synchro- 
nization. 

We also considered a periodic modulation of the pa- 
rameters, where the Eqn. 0can be replaced by 
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(4) 



FIG. 1: Synchronization is maintained in the presence of pa- 
rameter fluctuations. (j> — 1000 and Ap = 0.05 



where / is the frequency and a is the amplitude of 
modulation. By choosing an appropriate value for a and 
by changing / the quality of synchronization for various 
modulating frequencies can be studied. 



III. NUMERICAL EXAMPLES 

Coupled Rossler oscillators are a well known for nu- 
merical studies in synchronization. Consider a system of 
bidirectionally coupled Rossler oscillators. The coupled 
equation can be written as. 

x\ = -yi - z\ + c(x 2 - xi) (5) 

y\ = xi + pm 

z\ = 0.2 + zi(a;i - 10) 

X 2 = ~yi - Z2 + c{xi - X 2 ) 

2/2 = x 2 + p 2 y 2 

z 2 = 0.2 + z 2 (x 2 - 10) 



Here, the coupling strength c = 0.15, and po = 0.18. 
Though the coupling strength can also affect synchro- 
nization, we chose a value that is best suited for illus- 
trating the concepts. Also the value of Ap was fixed to 
be 0.05 for all fluctuation rates. 

FigDshows the synchronization plot in the presence of 
parameter fluctuations. It can be seen that the synchro- 
nization is robust. Fig[21shows that the coupled systems 
posses a parameter mismatch at any instant of evolu- 
tion of the system. Also at times the instantaneous mis- 
matches can be compared to the value of the average 
value parameter value itself. With the same value of Ap 
the synchronization is destroyed with a lower fluctuation 
rate as shown in FigOl 

We have studied the relationship between the syn- 
chronization error and the fluctuation rate. To quantify 
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FIG. 2: Synchronized system also posses a zero mean param- 
eter mismatch. 



the synchronization error we used the similarity function 
given by 



S 2 (r) 



{[x{t + r)-x{t)f) 

{( X 2(t))(x2(t))}h ■ 



(6) 



Here r is set to zero, which gives S(0), the error in 
synchrony. Fig. 0] Shows the plot of 5(0) vs. fluctuation 
rate. It can be seen the error diminishes rapidly with the 
increase in the fluctuation rate. 

A similar behavior was also seen with the periodic 
modulation of parameters. Fig. [S] shows that the syn- 
chronization error decreases with increase in the modu- 
lating frequency. This clearly suggests that the fluctua- 
tion frequency or rate is more important than the nature 
of fluctuations. 
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FIG. 3: Synchronization is destroyed in the presence of pa- 
rameter fluctuations^ with low fluctuation rates. Fluctuation 
rate (j) = 25 and Ap = 0.05 
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FIG. 5: The synchronization error decreases with the increase 
the modulating frequency. 
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FIG. 4: The synchronization error decreases with the increase 
in the fluctuation rate. 



IV. DISCUSSION 

The robustness of synchronization with high fluctu- 
ation rates and destruction of synchronization with low 
fluctuation rates can be explained analytically as follows. 
With Eqn^we can write an equation for the rate of sep- 
aration Xi — Xi of the trajectories as, 



d{Xi - X 2 
dt 



X 1 -X 2 =M( Pl ,p 2 ,X 1 ,X 2 ), (7) 



M(pi,p2, X±, X 2 ) is a function of the dynamical variables 
and the parameters of the coupled systems. This can be 
written as the sum of two terms, 



M(X 1 ,X 2 ) = S(X 1 ,X 2 ) + E(X 1 ,X 2 ). 



(8) 



ers or if the higher derivatives of SI w.r.t p is zero, 



Sl{p + x ) = fi (P) x ) + A P 



dSl(p,X) 
dp 



(9) 



This can be valid for functions in the dynamical equa- 
tions of a chaotic system if the parameter values are 
not near a bifurcation point. Using this relation, with 
Pi = Po + Api and p 2 = p + Ap 2 , 



S(p ,X 1 ,X 2 



and 



F 1 (p ,X 1 ,X 2 )-F i (po,X 1 ,X a ) 



E( Pl ,p 2 ,X 1 ,X 2 ) = A Pl — | P1=P0 



-Ap- 



dp 

dF 2 jp 1 , P2 ,X 1 ,X 2 ) 
dp 



with Api = £u and Ap 2 = £«. 

Here S(po,Xi,X 2 ) represents the quantity which of- 
fers a stable synchronization manifold, that is, when 
S(pi,p 2 , Xi,X 2 ) alone in the right hand side of the sepa- 
ration equation, coupled systems synchronize as t — > 00. 
The conditions for such a synchronization is widely dis- 
cussed in literature jl5| . The term E{pi 1 p 2l X\, X2) rep- 
resents the effect of the parameter mismatch. Coupled 
systems can synchronize if the overall effect of this term 
is zero as t — ► 00. One possible way for this is when 
E(pi,p 2 ,Xi, X 2 ,C) is of the form 



Eip 1 ,p 2 ,X 1 ,X 2 ,C) = 5^p(t) i i i (t) 



(10) 



This comes from the fact that for a function Slip, X), 
we can write for small Ap and neglecting its higher pow- 



where is the the fluctuation term and x(t) is 

the phase space variables of the coupled system. Here 
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the equation vanishes because the p(t)'s are zero mean 
rapidly fluctuating quantities and x(t)'s are the phase 
space variables that evolve slowly when compared to the 
the rapid fluctuations of the parameter. Thus x(t)'s 
can be assumed to be constant, in the time required 
for the fluctuations get summed to zero. This also ex- 
plains why the synchronization is destroyed when the 
fluctuation rate is low. With a low fluctuation rate 
the E(pi,p2, X\, X2, C) cannot be summed to zero ev- 
ery time since the phase space evolution time is compa- 
rable to the interval where a fixed parameter mismatch 
persists. Thus with a lower fluctuation rate the system 
always get time to respond to the parameter mismatch 
before it being canceled out. Also with a slowly varying 
parameter mismatch, a definite state of phase synchrony 
is also not attained and the system remains in a transient 
state through out the evolution in the phase space. 

In the present example, it can be seen that the quantity 
E(Pi,P2,Xi, X 2 ,C) can be expressed as, 

E(p 1 ,p 2 ,X 1 ,X 2 ,C) = & t |/i-6ti/a (11) 

because the fluctuating terms appear only in the equation 
of y only. Similar studies were done to the other param- 
eters of the coupled systems as well, which gave simi- 
lar result. Apart from gaussian random fluctuations, we 
studied perturbations with a uniform distribution. The 
results were qualitatively the same as for the gaussian 
perturbations which suggests that the most important 
quantity that determines the stability of synchronization 
is the fluctuation rates. 



V. CONCLUSIONS 

In this paper we studied the effect of parameter fluctu- 
ations on the synchronization of coupled chaotic systems. 
We investigated random parameter fluctuation and also a 



periodic modulation to the parameter. It was found that 
the most significant entity that determines the quality of 
synchronization is the fluctuation rates that we have de- 
fined or the frequency of fluctuation. Our study also show 
that the timescales with which the parameter fluctuates 
is more significant than the statistical or mathematical 
features of the fluctuations. 

The effect of noise on synchronization has been stud- 
ied in the past. Noise affects synchronization in different 
manner in various situations. In most of the cases noise 
destroys synchronization or make it unfit for the secure 
communication purposes that we have cited in the intro- 
duction 0, ^} . Also there are cases where synchroniza- 
tion is robust to noise [T^ or even induce synchronization 

The effect of noise and parameter fluctuations are dif- 
ferent. Noise induces perturbations to the phase space 
variables that decay while the system evolves. In a case 
where the parameter fluctuates, the resultant perturba- 
tions do not die out with the evolution of the system. 
It remains the same until it is corrected manually or 
the fluctuation modifies the parameter to a new value. 
Due to this reason, the fluctuation rate plays an im- 
portant role in determining the stability of synchroniza- 
tion in coupled chaotic systems. Parameter fluctuations 
may also have much higher significance in coupled arrays 
of nonlinear oscillators, especially in biological systems 
which exhibit synchronized behavior. Though we have 
not included these in our present paper, we hope that 
our studies will be a motivation in this direction. 
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